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Calculation  of  the  Derivative  of  the 
Parabolic  Cylinder  Function 


1  Introduction 

This  report  describes  the  numerical  evaluation  of  the  derivative  of  the  parabolic  cylinder  func¬ 
tion,  Dp(z),  where  Dp{z)  is  the  parabolic  cyUnder  function  defined  by  Whittaker  and  Watson  [1, 
pp.347-349].  (See  also  [2,  Ch.l9],  [3,  Ch.46],  and  [4].)  The  numerical  evaluation  is  for  complex 
argument  z  =  and  complex  order  p  =  — 1/2  +  iy,  x  and  y  real.  This  case  is  needed 

to  calculate  the  currents  excited  on  a  perfectly  conducting  parabolic  cylinder  by  a  transverse 
electric  (TE)  plane  wave,  and  the  corresponding  scattered  far  field  pattern  [5,  Eqs.  (7.47)  and 
(7.48)].  This  is  the  motivation  for  this  report,  designed  to  complement  the  valuable  communi¬ 
cation  of  Blanchard  and  Newman  [6]  which  treats  the  numerical  evaluation  of  the  same  case  of 
the  parabolic  cylinder  function  itself,  needed  for  the  corresponding  transverse  magnetic  (TM) 
scattering  calculations  [5,  Eqs.(7.12)  and  (7.13)]. 

2  Numerical  Evaluation  of  D^^z) 

We  present  expansions  for  Dp{z)  for  z  =  xa,  p  =  — 1/2  +  iy,  x  and  y  real,  and  a  =  The 

discussion  closely  parallels  that  of  Blanchard  and  Newman  [6].  We  adopt  their  division  of  the  xy 
plane  into  five  regions,  denoted  Regions  A-E,  as  shown  in  Figure  1,  and  give  numerically  useful 
expressions  for  calculating  Dp{z)  for  each  of  these  regions. 

Region  A:  |x|  <  5,  \y\  <  4. 

In  this  region,  the  exact  power  series  for  Dp{z)  can  be  used,  obtained  by  differentiating  the 
Received  for  publication  23  May  1996 
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exact  power  series  of  Dp{z)  given  in  [6]: 


D'{z)  = 


2-(p+5)/2g-z2/4 

tFS 


p 


A-p{-s/2zy 


+  i:  [2*^  -  2(f-  +  1)1  iV2zr 


r  (-2p  -  l){V2z)  +  ^ 

''  L  fe=i 


(2fc)! 

(-f)---(fc-l-i) 

(2fc  +  l)! 


[2k  +  l-2{p  +  l)]  (V2z)2'=+' 


(1) 


for  p  7^  0, 1,2, ...,  and  where  F  is  the  gamma  function  [2,  Ch.6].  On  a  computer  with  16  digit 
double  precision,  this  expression  gives  at  least  6  digit  accuracy  in  Region  A  if  the  summations 
are  continued  until  the  terms  decrease  below  1  •  10"^°  in  magnitude.  The  gamma  function  of 
complex  argument  is  also  computed  in  double  precision,  using  the  relations  (2)  and  (3),  derived 
using  the  recurrence,  reflection,  and  duplication  formulas  for  the  gamma  function  [2,  Ch.6],  to 
obtain  r(— p)  and  r(l/2  —  p/2)  from  r(— p/2); 


r(-;^) 


^(-1) 


COS  (f )  r-  (-1) 


(2) 

(3) 


where  *  denotes  complex  conjugation. 

Although  the  exact  power  series  (1)  for  Dp{z)  is  absolutely  convergent  for  all  x  and  y,  the 
terms  of  the  sums  grow  very  large  before  they  become  small  if  |a:|  or  |p|  is  large,  just  as  happens  for 
Dp{z)  itself  [6].  The  accuracy  of  (1)  is  then  degraded  by  the  process  of  subtracting  large  numbers 
of  comparable  magnitude,  and  also  eventually  by  the  generation  of  numbers  that  overflow  the 
dynamic  range  of  the  computer. 


Region  B:  x  >  5,  |p|  <  2. 

In  this  region,  the  asymptotic  series 
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Figure  1;  Definition  of  Regions  A-E  [6]. 


p{p  +  S)  p{p  -  l)ip  -  2){p  +  5) 

""  2  f  l!(2z2)i  2!(2z2)2 

P{P  -  1)(P  -  2)(p  -  3)(p  -  4)(p  +  7)  1  u\ 

3!(2z2)3  J  ^  ^ 

is  used.  This  series  can  be  obtained  by  differentiating  the  asymptotic  expansion  for  Dp{z)  given 
in  [1,  p.347]  or,  more  easily,  by  using  the  recursion  formula  [2,  Eq.  19.6.3]; 

i>;w  =  ]ipi>p-iW-.DpH.,(2)].  (5) 

The  expansion  (4)  is  valid  for  all p  and  for  |  arg(z)|  <  37r/4  whenp  <  z  [1,  p.347].  The  summation 
in  (4)  is  continued  until  the  terms  start  to  increase  in  magnitude.  At  least  five  significant  figure 
accuracy  can  be  obtained. 

Region  C:  All  x,  and  y>2  except  in  regions  A  and  E. 
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(6) 


In  this  region  we  use  a  large  order  asymptotic  expansion  obtained  by  Olver  [7]:  ^ 


+iy 


{xa) 


s=0 


i^Vs{t) 


where 

H'  =  ix'  >  0, 

t'  =  - ^ 

R  =  (t'2  +  1)1/2^ 


g{w)  =  2' 


-(l+n;2)/4p-iu2/4^-(l-w2)/2 


i  +  E 


92s+\ 
4s+2 


1-1 


W 


9l  =  53 


2021 

‘207360’ 


95  — 


s=Q 

5149591 

418037760’ 


^  +  ln(|t'|  +  R)] ,  sgn(i')  =  ||j, 


Vo  =  1,  5i  = 


r  -  et' 


15t'^  +  327t'2  -  143 


U2  = 


>2 


Vs 


24  ’  1152 

4042t'®  +  18189t"  +  36387t'®  +  238425t'2  -  259290t' 

414720  ■ 


The  terms  given  here  give  at  least  four  significant  figure  accuracy,  with  the  accuracy  improving 
as  y  increases.  Additional  terms  in  the  expansion  can  be  obtained  using  algorithms  given  by 
Olver  [7]. 


Region  D:  All  x,  and  y  <  -2  except  Regions  A  and  E. 


In  this  region  we 


use  the  following  large  order  asymptotic  expansion  obtained  by  Olver  [7]:  ^ 

‘bJ  ''  $  \  A'*' 

oo  n  "1 

-  >  0  (7) 

s=o  9-  ) 


^Equation  (6)  is  obtained  starting  with  [7,  Eq.(4.13)],  identifying  y  with  t  with  -ix/{2^),  and 

using  [7,  Eqs.(3.4),  (11.6),  (11.7),  (11.8),  (6.1),  (4.12),  (4.14),  (4.16)  ,  and  Fig.  (2b)].  See  also  [6,  Eq.(6)]  and  [8, 
Sect.  12.6.5). 

^Equation  (7)  is  obtained  starting  with  [7,  Eq.(8.15)],  identifying  y  with  ,  t  with  a:/(2y|^),  using 

[7,  Eqs.  (8.16), (7.4),  (11.21),  (8.18),  (8.19),  (4.14),  (4.16),  (8.13)],  and  [2,  Eq.  (10.4.9).  Equation  (8)  is  similarly 
obtained  starting  again  with  [7,  Eq.(8.15)],  letting  y  y^'^  ,  and  using  [7,  Eq.(5.24)].  See  also  [6,  Eqs.(7)  and 
(8)].  For  general  background  on  these  expansions  see  [9,  Ch.ll). 
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(8) 


D'.i+iy{xa) 


.C^ 


s=0 


IX'- 


where  Ai,  Bi,  Ai',  and  Bi'  are  the  Airy  functions  and  their  derivatives  [2,  Sect.  10.4],  and  gifx'a) 
is  defined  as  for  Region  C.  The  other  quantities  are  defined  by 


\/^> 


t  = 


kr 


V2/.i'’ 


=  t{e  - 1)'/2  -  iii[t + - 1)'/%  t  >  1,  (c  >  0), 
3 

1  (_C)3/2  =  arccos(t)  -  t{e  -  lf/\  -1  <  t  <  1,  (C  <  0), 


1/4 


$  = 


,t2-l 


2s+l 


'(^2  _  ^'^V2s-m+l{t), 


Do  =  l,  Ds  = 


1 


2s 

E 

7n=0 


■  _  l)3s  E  ®’"^2s-mW, 


,,  +  15t"  -  327t^  - 143 

”  1?  ^i(^)  24  ’  ^2w  1152  ’ 

-4042t®  +  18189t'^  -  36387t^  +  238425t^  +  259290t 


vsit)  = 


do  —  1) 


414720 

(2m  +  l)(2m  +  3)  •  •  •  (6m  —  1) 
m!(144)"* 


6m  + 1 

bm  =  ^  7  7 

6m  —  1 


For  |t  -  1|  <  0.25,  $  and  C  should  be  computed  by  the  expressions  given  in  [6],  while  for 
|t  -  1]  <  0.05,  the  following  expansions,  obtained  using  Mathematica,  should  be  used  for  Co,  Ci, 

and  Di- 

6  =  t-l 

$2 


Co  = 


j.  ^  r' 

-  +  ^S  + 


t  +  lV5  25 


97 

2250^ 


6^- 


173 

385000 


$2  /  88579 


1480656484897 


^  -i-  ,  88579  11626038701  2  i^ouuou-io^o^/ ^3  , 

=  (t  +  1)4  1,462000  220000*^  22702680000  4824319500000 


■)' 


5 


D, 


1 

{t+ir 


p041  48577  279481  ^2  21114659  3 

V  25200  ^  2772000*^  ^  840840000*^  ~  227026800000^ 


The  terms  given  here  give  at  least  four  significant  figure  accuracy  for  Region  D. 


Region  E:  x  <  —5,  —3.5  <y<  3.5. 


In  this  region  we  use  the  relation,  valid  for  all  z  and  p,  obtained  from  [1,  Sect.  16.52]: 
d;(z)  = 

For  p  =  —1/2  +  iy  and  z  = 


(9) 


(10) 


where  D'J—z)  is  calculated  as  above  for  Regions  B,  C,  or  D. 


3  Numerical  Results 

To  illustrate  the  accuracy  of  the  expansions  we  have  presented,  in  Table  1  we  compare  the 
results  obtained  with  the  expansions  for  Regions  B-E  with  the  Region  A  results  for  12  points 
on  the  perimeter  of  Region  A  (see  Figure  1),  corresponding  to  the  results  given  by  Blanchard 
and  Newman  [6]  for  the  parabolic  cylinder  function  itself.  Table  1  gives  the  values  of  x  and  y, 
the  region  whose  expansion  was  used  in  the  calculation,  the  calculated  value  of  Dp{z),  and  the 
number  of  significant  digits.  The  exact  power  series  results  for  Region  A  have  been  checked  with 
Mathematica  and  are  correct  to  the  number  of  figures  shown. 


4  Conclusion 

In  this  report  we  have  presented  numerically  useful  expressions  for  calculating  the  derivative  of 
the  parabolic  cylinder  function,  Dp{z),  for  the  special  case  p=  — 1/2  + iy  and  ;i:  =  This 

case  is  needed  for  calculating  the  currents  excited  on  a  perfectly  conducting  parabolic  cyhnder 
by  a  TE  plane  wave,  and  for  calculating  the  scattered  far  field  pattern. 
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Table  1:  Values  of  D'j,{z)  at  selected  points  in  Figure  1:  z  =  xe  and  p  =  -1/2  +  iy 


Point 

X 

y 

Region 

b;(z) 

Significant  Digits 

1 

0 

4 

A 

-0.214019E-f02  -  i0.867551E+01 

6 

1 

0 

4 

C 

-0.214038E+02  -  i0.867624E+01 

4.1 

2 

5 

4 

A 

-0.292821E+02  -  i0.58685lE-h00 

6 

2 

5 

4 

C 

-0.292820E+02  -  i0.586872E+00 

4.7 

3 

5 

2 

A 

0.557884E+01  -  i0.146969E+01 

6 

3" 

5 

2 

B 

0.557885E+01  -  i0.146966E-f-01 

5.5 

3 

5 

2 

C 

0.557886E+01  -  i0.146962E-h01 

5.2 

4 

5 

0 

A 

-0.103126E+01  +  i0.437434E+00 

6 

4 

5 

0 

B 

-0.103126E+01  -f-  i0.437433E-l-00 

6 

5 

5 

-2 

A 

0.197734E-h00  -  i0.775546E-01 

6 

5 

5 

-2 

B 

0.197733E+00  -  i0.775537E-01 

5 

5 

5 

-2 

D 

0.197735E+00  -  i0.775577E-01 

4.5 

6 

5 

-4 

A 

-0.384484E-01  +  i0.378112E-02 

6 

6 

5 

-4 

D 

-0.384484E-01  -I-  i0.378114E^02 

5.7 

7 

0 

-4 

A 

-0.214019E+02  -f  i0.86755lE-l-01 

6 

7 

0 

-4 

D 

-0.214033E-h02  +  i0.867600E-l-01 

4.3 

8 

-5 

-4 

A 

0.581438E+04  -  i0.235691E-F04 

6 

8 

-5 

-4 

D 

0.581437E+04  -  i0.235691E-h04 

6 

9 

-5 

-3.5 

A 

0.424382E+04  -  i0.236275E-l-03 

6 

9 

-5 

-3.5 

D 

0.424382E+04  -  i0.236275E+03 

6 

9 

-5 

-3.5 

E 

0.424382E+04  -  i0.236275E+03 

6 

10 

-5 

0 

A 

-0.190613E+01  -  i0.437434E+00 

6 

10 

-5 

0 

E 

-0.190613E-I-01  -  i0.437433E+00 

6 

11 

-5 

3.5 

A 

-0.112470E-I-02  -  i0.159690E+02 

6 

11 

-5 

3.5 

C 

-0.112468E-I-02  -  i0.159688E+02 

5.7 

11 

-5 

3.5 

E 

-0.112470E-I-02  -  i0.159690E+02 

6 

12 

-5 

4 

A 

-0.214258E-I-02  -  i0.199681E+02 

6 

12 

-5 

4 

C 

-0.214258E-I-02  -  i0.199680E+02 

6 
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MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab; 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 

The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications.  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


